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PREFACE 


This  report  was  prepared  by  Satoshi  Akagawa,  researcher  for  the  Shimizu  Corpora¬ 
tion,  Tokyo,  Japan,  while  he  was  a  visiting  scientist  with  the  Civil  and  Geotechnical  En¬ 
gineering  Research  Branch,  Experimental  Engineering  Division,  U.S.  Army  Cold  Regions 
Research  and  Engineering  Laboratory.  The  data  shown  in  this  report  were  observed  by 
the  author  at  the  Research  Institute  of  Shimizu  Corporation  in  1984  and  1985.  Most  of  the 
figures  in  this  report  were  acquired  from  the  data  analysis  performed  during  the  author's 
stay  at  CRREL  from  1986  to  1988.  During  the  analysis,  the  author  had  many  worthwhile 
discussions  with  the  CRREL  staff,  especially  Dr.  Richard  Berg,  Edwin  Chamberlain, 
William  Quinn,  Dr.  Yoshisuke  Nakano  and  Dr.  Patrick  Black.  The  author  also  appreciates 
the  assistance  of  David  Cate  in  helping  to  revise  this  report. 
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X-Ray  Photography  Method  for 
Experimental  Studies  of  the 
Frozen  Fringe  Characteristics  of  Freezing  Soil 

SATOSHI  AKAGAWA 


INTRODUCTION 

By  the  early  1980s,  the  study  of  frost  heave  appeared  to  be  complete  because  several 
numerical  frost  heave  models  had  been  developed.  Some  of  them  (Gilpin  1980,  O'Neill 
and  Miller  1982)  predict  realistic  ice-banding  patterns.  However,  these  models  are  based 
on  assumptions  that  still  remain  to  be  proven  by  precise  experiments.  Most  of  the 
assumptions  that  affect  the  reliability  of  the  models  are  related  to  the  "frozen  fringe," 
because  phenomena  that  occur  in  this  area  seem  to  initiate  the  entire  frost  heaving  process. 
Unfortunately,  insufficient  knowledge  about  the  frozen  fringe  makes  it  difficult  to  verify 
the  frost  heave  models. 

The  main  objectives  of  this  report  are  ' 

•  To  describe  the  development  of  an  X-ray  photographic  method  for  observing  the 
frost-heaving  properties  of  soil,  especially  the  characteristics  of  the  frozen  fringe; 

•  To  describe  observations  of  the  temperature  and  strain  field  change  in  the  frozen 
fringe  using  the  procedure  developed;  and 

•  To  publish  the  test  results  as  a  data  base  that  can  be  used  in  evaluating  frcst  heave 
models. 

The  method  used  X-ray  photography  with  thermocouples  installed  inside  lead  mark¬ 
ing  spheres  to  observe  soil  expansion,  including  expansion  due  to  in-situ  freezing  and 
expansion  due  to  segregation  freezing.  The  method  is  an  improvement  over  the  methods 
of  Penner  (1975),  Ishizaki  (1985)  and  Svec  (1986),  who  used  X-ray  radiography.  The  im¬ 
provements  include  an  image-processing  technique  for  analyzing  the  coordinates  of  the 
lead  spheres  on  the  photos,  a  three-dimensional  analysis  of  the  lead  spheres'  coordinates 
to  compensate  for  the  geometrical  distortion  caused  by  the  use  of  a  point  X-ray  source,  and 
direct  measurements  of  the  temperature  at  each  lead  sphere  using  thermocouples  in  the 
spheres.  Therefore,  both  the  deformation  and  the  temperature  field  in  a  freezing  soil  can 
be  observed  directly.  . 

The  experiments  produced  results  on  the  following  characteristics: 

•  Temperature  profiles  in  two  freezing  soil  specimens; 

•  The  expansion  activity  (strain  and  strain  rate)  distributions  in  two  specimens; 

•  The  apparent  temperature  range  of  the  frozen  fringe; 

•  The  apparent  thickness  of  the  frozen  fringe;  and 

•  The  thermal  conductivity  of  the  frozen  fringe  while  the  warmest  ice  lens  is 
segregating. 


1.  Specimen 

2.  Warm  plaie 

3.  Cold  plaie 

4.  Porous  plaie 

5.  Double  acrylic  cylinder 

6.  Dead  weight 


7.  Cooling  unit 

8.  Brine  temperature  controller 
0.  Burette  for  water  table 

10.  Deformation  transducer 

11.  Thermocouples 
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These  data  can  be  applied  to  studies  of  the  heaving  activity  distribution  in  conjunction 
with  temperature  field,  the  dominance  of  the  heave  mode  (in-situ  freezing  and  segrega¬ 
tion  freezing)  (discussed  in  Akagawa  1988a),  and  water  sources  for  segregating  ice 
(discussed  in  Akagawa  1988a).  The  method  could  also  be  used  to  study  the  permeability 
of  the  frozen  fringe.  The  present  technology  is  not  adequate  for  microscopic  observations, 
such  as  the  distribution  of  ice  and  unfrozen  water  around  soil  skeleton,  the  mechanical 
properties  of  absorbed  water,  and  the  ice  segregation  mechanism  in  freezing  soil. 


EXPERIMENTAL  METHODS 

Two  step-temperature-change  frost 
heave  tests  were  conducted  (Test  A  and 
Test  B  in  Table  1).  During  the  tests,  136  X- 
ray  photos  were  taken.  The  X-ray  photos 
were  taken  in  pairs  with  different  locations 
of  the  X-ray  source  so  that  the  three-di¬ 
mensional  coordinates  of  the  lead  spheres 
could  be  analyzed. 

Test  apparatus 

The  frost  heave  test  apparatus  is  shown  in  Figure  1.  It  has  one-dimensional  heat  flow 
through  the  cylindrical  specimen,  accomplished  by  reducing  the  heat  flow  normal  to  the 
specimen  axis  using  an  annular  dead  air  space  maintained  at  0°C.  The  overburden 
pressure  was  kept  constant  by  using  a  dead  weight  for  loading.  The  apparatus  has  a 
ripple-free  temperature  control  employing  a  large-volume  (11  L)  circulating  coolant 
system.  It  uses  drift-free  deformation  transducers  for  recording  the  total  heave,  accom¬ 
plished  by  using  a  digital  deformation  measuring  system  without  an  amplifier.  The  rest 
of  the  features  of  this  frost  heave  cell  are  given  by  Akagawa  (1983, 1985). 

The  X-ray  facilities  (Fig.  2)  have  the  appropriate  X-ray  radiation  strength  obtained  by 
applying  100  kV  and  5  mA  to  the  X-ray  tube.  The  X-ray  source  can  move  between  two 


Table  1.  Test  conditions. 


Test  A 

TesIB 

Overburden  pressure  (kPa) 

60 

no 

Pore  water  pressure  (kPa) 

0 

0 

Cold-side  temperature  (°C) 

-5.9 

-5.5 

Warm-side  temperature  (°C) 

2.2 

3.0 
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positions  in  a  plane  parallel  to  the  X-ray  film,  as  shown  in  Figure  3.  This  enables  X-ray 
photos  to  be  taken  from  two  consistent  (but  different)  X-ray  source  positions,  allowing  a 
three-dimensional  coordinate  analysis.  The  difference  in  X-ray  path  length  through  the 
circular  section  of  the  specimen  is  compensated  for  by  using  the  axisymmetrical  alumi¬ 
num  wedge  filter  shown  in  Figure  2  (this  technique  was  used  only  in  Test  B). 

The  lead  spheres  and  thermocouples  are  shown  in  Figure  4.  The  timing  and  position 
of  the  coordinate  and  temperature  measurements  matched  precisely  because  the  thermo- 
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Height  from  the  Cold  Plate  in  Test  B  (mm) 


couples  were  installed  within  the  lead  spheres.  The  thermal  disturbance  created  by  the 
installation  of  the  spheres  and  thermocouples  was  negligible  because  the  spheres  were 
very  small  (1.5  mm  in  diameter)  and  the  sensors  were  thin  (0.01  mm  in  diameter). 

The  thermocouples  were  calibrated  by  tying  them  into  a  bundle  and  putting  them  into 
amix  tureofcooled  water  and  crushed  iceatatemperatureofaboi’<--200C.Thebundlewas 
left  until  the  temperature  of  the  mixture  rose  above  the  freezing  point.  Then  the  recorded 
temperatures  of  the  thermocouples  at  which  water  and  ice  were  in  equilibrium  were  read, 
and  the  thermocouples  that  showed  0.0  C  during  the  two-phase  equilibrium  were 
selected .  The  temperatures  were  measured  with  a  digitally  linearized  data  logger  that  had 
been  calibrated  at  the  factory  to  an  accuracy  of  0.1°C.  All  the  thermocouples  worked 
during  the  frost  heave  test  except  Numbers  7  and  8  in  Test  A.  Thermocouple  Number  7 
stopped  functioning  after  about  200  hours  and  Number  8  stopped  after  500  hours. 


Soil  specimens 

Some  physical,  mechanical  and 
thermal  properties  of  the  soil  stud¬ 
ied  are  listed  in  Table  2.  The  soil 
was  overconsolidated  and  con¬ 
tained  a  lot  of  unfrozen  water.  This 
test  specimen  was  cored  from  a 
water-saturated  alluvium  bed  and 
trimmed  in  the  undisturbed  state. 
A  real  undisturbed  soil  was  used  so 
that  real  frost  heave  phenomena 
would  be  observed.  The  result  and 
applicability  of  proposed  frost 
heave  models  can  be  verified  by 
comparing  them  with  the  experi¬ 
mental  data  shown  in  this  report. 

Test  conditions 

The  test  conditions  are  listed  in 
Table  1.  The  temperatures  at  both 
the  top  and  the  bottom  surfaces 
were  kept  constant  exceptbetween 
360  and  377  hours  and  930.5  and 
1006  hours  in  Test  A  and  between 
111  and  167  hours  in  Test  B.  The 
overburden  pressures  were  also 
kept  constant  up  to  329  hours  in 
Test  A  and  383  hours  in  Test  B. 
During  the  growth  of  the  final  ice 
lens  (stationary' heaving),  the  warm- 
side  temperature  T.  the  cold-side 
temperature  T  and  the  overbur¬ 
den  pressure  P  were  changed  sev¬ 
eral  times  to  observe  the  effects  on 
heave  rate. 


Table  2.  Physical,  mechanical  and  thermal  prop 
erties  of  the  soil. 


Size  distribution  (%  by  weight) 


Clay  (<0.002  mm) 

45.0 

Silt  (0.002-0.074  mm) 

41.0 

Sand  (0.074-2.000  mm) 

14.0 

Gravel  (>2.000  mm) 

0 

Atterberg  limits 

LL 

101.0(%) 

IP 

68.7 

Specific  gravity 

2.48 

Wet  density 

1.45  g/cm3 

Void  ratio 

2.6 

Preconsolidation  stress 

0.248  MPa 

Consolidationcoefficient 

5.78xl0~2cm2/min 

Compression  index 

1.233 

Permeability 

By  consolidation  test 

1.77x10  cm/min 

By  permeability  test 

4.61  x  1C6  cm/min 

Volumetric  unfrozen  water  content  by  TDR 

at  -0.2°C 

60.7% 

-0.3 

53.2 

-0.4 

38.0 

-0.6 

29.8 

-0.75 

27.7 

-1.0 

25.1 

-1.5 

22.1 

-2.0 

20.2 

-3.0 

18.3 

-5.0 

14.6 

-7.0 

13.1 

-10.0 

12.3 

Thermal  conductivity  by  thermal  probe  method 

at  24°C 

6.0cal/cm-hr-°C 

-5 

12.0 

-10 

16.1 

In  this  report  the  frost  heave  modes  of  the  s  tep-temperature-change  frost  heave  test  are 
identified  as  "transient  heaving"  and  "stationary  heaving"  (Akagawa  1985).  Transient 
heaving  is  characterized  by  soil  freezing  that  causes  a  deu  ease  in  the  unfrozen  soil  length, 
i.e.  a  penetrating  0°C  isotherm.  It  results  in  what  is  called  ice  cement  first  and  then  ice 
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banding.  Stationary  heaving,  on  the  other  hand,  is  characterized  by  soil  freezing  that 
causes  no  decrease  in  the  unfrozen  soil  length,  i.e.  a  stationary  0°C  isotherm.  It  is  the  result 
of  the  growth  of  what  is  called  the  last  or  final  ice  lens. 

Test  procedure 

The  frost  heave  test  was  conducted  using  the  following  procedure.  The  undisturbed 
specimen  was  placed  in  the  water-filled,  double-acrylic  cylinder  on  top  of  the  bottom  cold 
plate.  (The  test  sample  was  frozen  from  the  bottom  upward.)  After  the  warm  plate  was 
placed  on  top  of  the  specimen,  the  air  bubbles  were  extracted  from  the  space  between  the 
top  of  the  soil  and  the  warm  plate  through  the  upper  pore-water  supply  line.  An 
overburden  pressure  of  100  kPa  was  applied  in  addition  to  the  test  pressure  to  ensure  good 
contact  between  the  temperature  control  plate  and  the  specimen  surface.  This  pressure 
was  maintained  until  consolidation  was  complete.  During  consolidation  the  warm-plate 
and  cold-plate  temperatures  were  maintained  at  +3°C.  The  overburden  pressure  was 
reduced  to  the  value  shown  in  Table  1  by  unloading  the  extra  100-kPa  surcharge  and 
waiting  until  the  rebound  was  complete.  A  cold  brine  with  enough  heat  removal  capacity 
to  cool  the  cold  plate  and  brine  circulating  system  to  the  specified  temperature  was 
circulated  to  start  the  test,  and  then  the  surface  temperature  of  the  cold  plate  was  adjusted 
to  the  test  temperature  when  the  temperature  at  the  cold  plate  dropped  slightly  below  the 
test  temperature. 

X-ray  photograph  exposure  procedure 

The  X-ray  photos  were  taken  at  the  times  listed  in  Table 
3.  For  the  first  exposure  the  X-ray  film  was  set  behind  the 
specimen  at  a  right  angle  to  the  X-ray  source  through  the 
center  of  the  specimen,  shown  as  Position  1  in  Figure  3.  The 
specimen  was  then  exposed  to  X  rays  for  3-5  minutes 
depending  upon  the  thickness  of  the  ice  lens.  Another  X-ray 
film  was  set  in  place,  and  the  X-ray  source  was  moved  from 
Position  1  to  Position  2,  shown  in  Figure  3.  A  second  expo¬ 
sure  was  taken  using  the  same  procedure  as  the  first. 

A  typical  X-ray  photo,  taken  at  1 18  hours  during  Test  B, 
is  shown  in  Figure  5  along  with  a  temperature  profile.  All 
the  X-ray  photos  are  included  in  Appendix  G.  Since  each 
lead  sphere  appears  clearly  in  the  photo,  the  coordinates  of 
each  sphere  can  be  measured  precisely.  Since  the  lead 
spheres'  temperatures  have  also  been  measured,  a  precise 
temperature  profile  associated  with  each  photo  is  available. 

The  ice  lens  distribution  is  also  visible  in  this  X-ray  photo. 

Therefore,  the  temperature  where  the  warmest  ice  lens  is 
growing  can  be  determined. 

The  X-ray  photo  in  Figure  5  shows  20  lead  spheres.  Fifteen  of  them,  numbered  0  to  14, 
represent  the  lead  spheres  that  were  installed  within  the  soil  specimen.  The  rest  of  the 
spheres,  numbered  from  21  to  25,  were  set  between  the  soil  specimen  and  the  X-ray  film; 
these  five  X-ray  shadows  were  used  as  fixed  control  points  for  the  coordinate  analysis  for 
Test  B.  For  Test  A,  four  immobile  points  numbered  21  to  24  were  used. 

Image  processing  for  the  X-ray  photographs 

An  image  processing  procedure  was  used  so  that  the  coordinates  of  the  lead  spheres 
could  be  determined  by  computer.  The  X-ray  photo  image  was  converted  to  640  x  400 
pixels  of  digital  data  using  a  high-resolution  CCD  camera  that  has  8-bit  (256)  intensity 
grades  (one  pixel  represents  an  area  of  0.17  x  0.17  mm  on  the  X-ray  photograph). 


Table  3.  Hours  when  X- 
ray  photographs  were 
taken. 


Test  A  Test  B 


-0.5  to  23.5* 

-1  to  24’ 

44.5 

28 

141.5 

46 

304.5 

70 

524.5 

118 

664.5 

166 

832.5 

172 

1408.5 

214 

286 

382 

388 

502 

646 

*  Every  hour  between  these 
times. 
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The  digitized  data  were  processed  by  two  kinds  of  window  operators  (spatial  filtering). 
First,  a  3  x  3  low-pass  filter  was  applied  to  improve  the  signal-to-noise  ratio,  and  then  a 
3x3  spatial  differential  filter  was  applied  to  enhance  the  edge  of  the  lead  spheres'  images. 
By  this  process  the  image  of  a  lead  sphere  was  convened  from  a  low-intensity  black 
circular  spot  to  the  ring  pattern  shown  in  Figure  6a.  A  computer  can  distinguish  this 
pattern  from  other  patterns,  such  as  vertical  shrinkage  cracks  and  lamina -chaped  ice 
lenses. 

The  addresses  or  coordinates  of  the  interior  of  the  circle  from  which  a  centroid  can  be 
calculated  were  determined  by  applying  two  other  filters  (Fig.  6b, c).  Because  of  the  large 
number  of  pixels  in  a  photo  (256,000),  rough  filtering  requiring  rather  short  computation 
times  was  applied  first.  This  filter  (Fig.  6b)  searches  the  area  for  the  ring-shaped  image 
surrounding  the  lead  spheres  by  checking  nine  pixel  intensities.  If  a  pixel  has  low  intensity 
and  is  surrounded  by  four  outer  or  inner  high-intensity  pixels,  this  pixel  is  selected  for 
further  filtering.  The  subsequent  filterings  need  much  more  computation  time  than  the 
first.  These  are  only  applied  to  the  pixels  nominated  by  the  first  simple  filter.  These  filters 
(Fig.  6c)  determine  the  shape  of  the  ring  by  checking  the  integrated  intensity  of  the  pixels 
in  the  high-intensity  ring  and  the  low-intensity  circle  in  the  center.  Only  about  five  pixels 
located  in  the  interior  of  the  ring  pass  all  these  filtering  stages.  The  centroids  of  all  the  ring 
images  were  computed  by  calculating  the  mean  value  of  every  nominated  pixel's 
coordinates. 

Since  these  coordinates  represent  only  relative  locations,  the  absolute  coordinates  were 
computed  by  comparing  them  to  a  reference  coordinate.  In  this  conversion  the  coordi¬ 
nates  of  lead  sphere  22  in  X-ray  photos  taken  from  Position  1  were  selected  as  the 
reference.  The  three-dimensional  coordinates  of  the  ring  image  were  calculated  using  the 
geometry  shown  in  Figure  3.  The  mean  value  of  the  coordinates'  standard  deviation  for 
the  immobile  points  (numbers  21, 23, 24  and  25)  was  found  to  be  0.053  mm. 
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c.  Second  filter. 


Figure  6.  Image-processing  filters  for  analyzing  the  sphere  coordinates. 
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Figure  7.  General  Iwave  prop¬ 
erties  of  the  soil-  A,  B  and  C 
are  constants,  and  dh/dt  is 
the  total  heave  rate.  (After 
Akagawa  1985.) 


TEST  RESULTS 
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Conventional  heave  test  results 

Total  heave  and  total  heave  rate  curves  of  this  soil  under  step-temperature-change 
freezing,  discussed  by  Akagawa  (1983. 1985),  are  shown  in  Figure  7.  As  in  other  similar 
frost  heave  tests,  this  soil  consistently  heaves  at  a  high  rate  during  the  early  stages  of  freez¬ 
ing  with  one  minimum  point  and  at  a  gradually  decreasing  rate  during  the  later  stages. 

Figure  8  shows  that  the  same  types  of  heaving  properties  were  observed  in  both  Test 
A  and  Test  B.  The  figure  also  shows  the  durations  of  the  overburden  pressure  and  the 


Time  (hr) 


a.  Test  A. 


X:  Uncertain  disturbance  of  the  heave 
rate,  possibly  caused  bi/  disturbance  of  the 
coldroom  temperature. 

A:  Heave  rate  increase  due  to  lowering  of 
the  warm-side  temperature. 

B:  Heave  rate  decrease  due  to  raising  of 
the  warm-side  temperature. 


C:  Heave  rate  increase  due  to  lowering 
of  the  cold-side  temperature. 

D:  Heave  rate  decrease  due  to  raising  of 
the  cold-side  temperature. 

E:  Heave  rate  decrease  due  to  raising  of 
the  overburden  pressure. 

F:  Heave  rate  increase  due  to  lowering  of 
the  overburden  pressure. 


Figure  8.  Conventional  results  of  frost  heave  tests. 
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b.  Test  B. 

Figure  8  (cont'd). 

times  when  X-ray  photos  were  taken.  In  Figure  8b  (Test  B)  the  water  intake  amount  is 
shown  as  the  height  of  water  in  the  cell  (the  amount  of  intake  water  divided  by  the 
horizontal  section  area  of  the  cell).  Figure  9  shows  the  water  intake  amount  for  Test  A. 
During  the  early  stages  of  Test  A  (13-14  hours),  pore  water  started  to  leak  from  a  small  hole 
in  the  upper  part  of  the  cell  through  which  a  thermocouple  was  installed.  After  that  the 
water  head  of  the  pore  water  source  was  maintained  just  below  the  level  of  the  hole  to 
prevent  leakage.  Thus  the  water  flow  measurements  for  Test  A  were  abandoned  after  13 
hours,  when  the  leak  began. 

Figure  8  confirms  the  following  experimental  results: 

•  Under  similar  thermal  conditions,  the  lower  overburden  pressures  resulted  in 
higher  total  heave  amounts.  (The  total  heave  amount  in  Test  A  was  greater  than  in 
Test  B.)  However,  this  trend  is  not  clear  during  transient  heaving. 

•  Raising  the  cold-plate  temperature  at  377  hours  in  Test  A  and  111  hours  in  Test  B 
resulted  in  a  decrease  in  heave  rate,  and  lowering  the  cold-plate  temperature  at  360 
hours  in  Test  A  and  167  hours  in  Test  B  increased  the  heave  rate. 


11 


Figure  9.  Pore  water  leakage  after  13  hours  in  Test  A. 


*  The  heave  rate  varied  during  step-temperature  freezing.  Therefore,  a  simple  de¬ 
scription  of  the  heave  rate  or  heave  amount  seems  meaningless. 

These  are  the  experimental  results  normally  observed  by  conventional  measurements 
during  frost  heave  tests. 

Positions  of  the  lead  spheres 

The  changes  in  the  Z  coordinates  of  all  the  lead  spheres  in  Test  A  and  Test  B  are  shown 
in  Figure  10  in  comparison  with  0°C  isotherm.  Each  data  point  represents  the  Z  coordinate 
or  the  distance  from  the  cold  plate  at  a  specific  time.  Each  distance  between  two  adjacent 
data  points  at  a  specific  time  shows  the  thickness  of  the  soil  layer  between  them.  This 
figure  shows  when  and  where  the  soil  expands.  As  the  0°C  isotherm  penetrates,  the  soil 
layers  below  it  expand.  However,  since  it  is  difficult  to  make  precise  observations  from 
these  figures,  the  expansion  properties  of  each  soil  layer  will  be  shown  later  in  detail. 

Since  the  0°C  isotherm  penetration  into  the  unfrozen  soil  seems  to  become  stationary 
at  about  40-50  hours,  stationary  heaving  had  presumably  started  around  this  time. 

Temperature  profiles 

From  the  Z  coordinates  and  the  temperatures  of  the  lead  spheres,  temperature  profiles 
can  be  drawn  (Fig.  11).  Appendix  A  lists  the  Z  coordinates  and  temperatures  at  each  ob¬ 
servation  time  for  both  tests.  All  the  temperature  profiles  are  shown  in  Appendices  C  and 
E  along  with  die  strain  or  strain  rate  distributions. 

Each  temperature  profile  shows  the  location  of  the  warmest  ice  lens  if  the  lens  was 
visible  in  the  X-ray  photo.  Since  past  definitions  of  the  term  "frozen  fringe"  are  not  clear, 
the  distance  between  the  location  ol  the  visible  warmest  ice  lens  and  the  0°C  isotherm  will 
be  assumed  to  be  the  apparent  thickness  of  the  frozen  fringe  Lff.  The  temperature  at  the 
warm  side  of  the  warmest  ice  lens  is  assumed  to  be  the  apparent  ice  segregation  tem¬ 
perature  T . 

During  Test  A,  temperatures  at  the  side  bo  andary  between  the  frost  heave  cell  and  the 
test  specimen  were  also  measured  every  5  mm.  Figure  12  shows  the  temperature  profiles 
for  the  side  of  the  specimen  and  the  center  of  the  specimen.  Even  though  care  was  taken 
to  prevent  it,  lateral  heat  flow  may  have  caused  a  slight  temperature  difference,  especially 
at  the  lower  end  of  the  cell.  This  tendency  is  reasonable  because  the  temperature  difference 
at  the  cold  end  of  the  specimen  is  the  largest;  the  room  temperature  was  kept  at  about  0°C 
and  the  specimen  temperature  was  kept  at  about  +2°C  at  the  warm  side  and  -6°C  at  the 
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Figure  11.  Analyzed  temperature  profiles. 


a.  Test  A. 


b.  Test  B. 
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a.  After  0.5  hours. 


b.  After  5.5  hours. 

Figure  12.  Comparison  of  temperature  profiles  for  the  lead  spheres  and  the  surface 
of  the  samples. 
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Temperature  l°C)  Temperature  (°C) 


c.  After  10.5  hours. 


d.  After  15.5  hours. 

Figure  12  (cont'd).  Comparison  of  temperature  profiles  for  the  lead  spheres  and  the 
surface  of  the  samples. 
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Figure  13.  Strain  development  in  freezing  soils.  Heave  is  positive. 
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Figure  13  (cont'd). 


cold  side.  Since  the  ice  lenses  in  the  specimen  were  not  concave  but  were  rather  flat  this 
lateral  heat  flow  was  not  causing  enough  thermal  disturbance  to  seriously  affect  ice  lens 
formation,  especially  at  the  warmest  ice  lens,  where  segregation  was  presumably  most 
active.  However,  these  comparisons  suggest  that  the  measurement  error  has  to  be  recog¬ 
nized  when  temperature  profiles  are  measured  from  the  side  surface  of  the  specimen. 

A  much  more  fundamental  problem  is  whether  or  not  linear  interpolations  of  the 
temperature  profile  from  the  warm-side  and  cold-side  temperatures  are  acceptable.  Of 
course  the  validity  of  the  interpolation  will  depend  on  how  precise  a  temperature  reading 
is  necessary  for  the  specific  study.  However,  from  the  temperature  profiles  in  Figure  12, 
the  linear  interpolation  seems  to  be  unacceptable  during  the  early  part  of  the  transient 
freezing  and  most  of  the  stationary  freezing  if  thick  ice  lens  segregation  is  expected  or 
recorded. 

Strain  developed  by  freezing  in  each  soil  layer 

Figure  13  shows  the  expansion  change  in  percent  strain  for  each  layer.  To  calculate 
strain,  two  adjacent  lead  spheres  at  time  f  were  identified  as  (t,n)  and  (t,n+l).  Tire  Z 
coordinates  for  these  spheres  were  written  as  Z(f,«)  and  Z(t,«+1).  The  conventional  strain 
e  in  a  soil  layer  was  then  calculated  as  follows: 


g  (n  t)  _  [ZM  -Z(f ,«+!)]  -[Z(0 ,«)  -  Z(0,?h-1)] 
Z(0,m)  -  Z(0,n+1) 


(1) 


Expansive  strain  is  assumed  to  be  positive.  Appendix  B  shows  the  strain  change  in  each 
soil  layer,  and  Appendix  C  shows  the  strain  distribution  at  times  when  X-ray  photos  were 
taken. 

In  Test  A,  since  the  strain  at  0.5  hours  is  treated  as  the  initial  value,  the  strain  developed 
between  0  and  0.5  hours  is  ignored;  0  hours  is  treated  as  the  initial  value  in  Test  B.  Because 
of  this  difference,  no  strain  was  shown  in  Layer  15  of  Test  A,  located  between  Sphere  14 
and  the  cold  plate,  because  the  temperature  of  Sphere  14  had  already  dropped  to  -2°C  at 
0.5  hours  and  about  70%  of  the  water  had  frozen  by  then.  Since  only  1-2%  of  the  unfrozen 
water  would  be  expected  to  freeze  during  a  further  temperature  drop  in  this  layer,  no 
expansion  was  measured  during  the  experiment. 

These  experimental  results  show  that  expansive  (heaving)  strain  progressed  from  the 
lower  layers  to  the  upper  layers  as  time  elapsed  following  penetration  of  the  0°C  isotherm. 
As  the  0°C  isotherm  penetration  rate  decreased,  expansive  strain  became  greater,  and 
negative  (consolidation)  strain  in  the  unfrozen  soil  appeared  clearly  before  the  expansion 
took  place. 

All  of  the  unfrozen  soil  layers  consistently  shrank  in  thickness  during  the  first  freezing 
stage.  This  kind  of  shrinkage  has  been  observed  by  Taber  (1929)  and  studied  by  Nixon  and 
Morgenstern  (1973)  and  Chamberlain  and  Gow  (1979);  however,  this  report  might  be  the 
first  description  of  when,  where  and  how  much  shrinkage  is  taking  place.  The  earlier 
studies  concluded  that  this  shrinkage  may  be  caused  by  the  increase  in  effective  pressure 
caused  by  the  decrease  in  pore  water  pressure  at  the  freezing  front.  The  magnitude  of  the 
shrinkage  observed  in  the  unfrozen  layers  in  this  study  tended  to  be  greater  in  the  layer 
nearest  the  warmest  ice  lens  than  in  the  layer  nearest  the  warm  plate.  This  tendency 
supports  the  conclusions  of  Taber,  Nixon  and  Chamberlain  that  there  is  a  pore  water 
depression  at  the  boundary  of  the  warm  side  of  the  warmest  ice  lens.  The  maximum 
shrinkage  is  so  high  (about  10%  in  strain)  that  it  results  in  much  greater  preconsolidation 
stress  of  the  clay.  The  preconsolidation  stress  for  this  soil  has  been  observed  at  a  strain  of 
about  5%  by  consolidation  tests  (Fig.  14). 
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Figure  14.  Consolidation  properties  of  the  test  soil  and  nomographic  pore  water 
pressure  of  segregating  ice. 

The  maximum  shrinkage  strain  of  10%  observed  in  this  study  would  be  caused  by  an 
increase  of  effective  pressure  of  about  0.4  MPa.  This  increase  of  effective  pressure  may  be 
understood  to  be  the  result  of  negative  pore  water  pressures  in  the  unfrozen  layers.  Since 
the  magnitude  of  the  pore  water  pressure  depression  is  the  mean  pore  pressure  in  the 
layer,  and  the  pore  water  pressure  may  have  a  gradient,  the  maximum  depression  should 
be  greater  than  the  revealed  value. 

In  addition,  even  though  the  shrinkage  had  been  studied  earlier,  many  frost  heave  tests 
have  been  conducted  without  concern  for  this  aspect.  In  other  words,  frost  heave  studies 
using  total  heave  amount  or  heave  rate  may  cause  confusion  because  the  total  heave 
measured  may  consist  of  both  heave  and  shrinkage,  particularly  in  clay  soils.  Figure  15 
shows  the  positive  expansion  (real  heave  amount)  and  the  negative  expansion  (real 
shrinkage)  observed  in  the  tests.  The  time-dependent  properties  of  the  real  heaving  and 
shrinkage  in  Test  A  and  Test  B  are  identical  qualitatively  and  quantitatively  if  we  consider 
the  initial  time  difference  mentioned  above.  The  strong  shrinkage  took  place  during  the 
first  3-4  hours,  and  a  slight  recovery  was  seen  between  4  and  20  hours.  This  drastic 
shrinkage  in  the  early  stages  of  freezing  is  clearly  reflected  in  a  negative  peak  in  the  total 
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b.  All  hours. 

Figure  25.  Expansion  (heave)  and  shrinkage  (consolidation)  in  the  freezing  soil. 
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heave  rate  curves  shown  in  Figure  16,  which  shows  the  heave  properties  in  the  early  part 
of  freezing.  This  negative  peak  has  been  observed  even  in  freezing  rocks  (Akagawa  et  al. 
1988). 

Therefore,  without  knowing  the  amount  of  shrinkage,  the  total  heave  amount  will 
demonstrate  nothing  for  this  time  period.  Of  course  this  time  period  and  the  magnitude 
of  shrinkage  may  vary  with  test  conditions  and  soil  properties.  This  period  is  related  to 
a  period  of  transient  freezing  that  is  characterized  by  a  penetration  of  the  0°C  isotherm. 
In  other  words,  the  total  heave  observed  during  stationary  frost  heaving  that  is  character¬ 
ized  by  the  dominant  growth  of  the  final  ice  lens  may  show  the  real  growth  amount  of  the 
ice  lenses,  as  is  seen  in  Figure  15b,  because  no  shrinkage  is  seen  during  this  period.  In 
addition,  the  effect  of  the  overburden  pressure  on  the  heaving  amount  becomes  clear 
during  stationary  frost  heaving,  but  the  effect  is  not  clear  during  transient  heaving  (Fig. 
15b). 

Strain  rate  of  each  soil  layer 

Figure  17  shows  the  time-dependent  expansion  rate  change  in  each  layer.  To  calculate 
the  strain  rate,  two  adjacent  lead  spheres  at  time  f  were  identified  as  (t^n)  and  (f r«+ 1 ).  The 
Z  coordinates  for  these  spheres  were  written  as  Z{tvn)  and  Z(/?/h+ 1 ).  The  mean  true  strain 
rate  at  a  period  between  t1  and  /,  (de/d/),  is  defined  as 


ds  ((/],«)  _  [Z  (/ 2, »)  -  Z  (t2,  n+ 1)]  -  [Z  (tv  h)  -  Z  (th  >h-1)]  (2) 

dt  [Z(#t, «)  -  Z  (*,,  »«-l)]  - 

Expansive  strain  rate  is  assumed  to  be  positive.  Appendix  D  shows  the  strain  rate  change 
in  each  soii  layer,  and  Appendix  E  shows  the  strain  rate  distribution  at  times  when  X-ray 
photos  were  taken. 

Figure  17,  a  smoothed  version  of  Appendix  D,  clearly  shows  the  deformation  activity 
(heave  and  consolidation)  due  to  freezing.  The  expansive  strain  began  in  the  lower 
(colder)  layers  at  the  early  stages  of  freezing  and  moved  to  the  upper  layers  as  the  freezing 
zone  reached  each  layer.  In  the  early  stages  of  freezing,  unfrozen  Layers  6-1 1  in  Test  A  and 
9-13  in  Test  B  shrunk  strongly.  Layers  1-5  in  Test  A  and  1-8  in  Test  B,  which  remained 
unfrozen  thro-  ghout  the  tests,  shrank  slightly  during  the  first  20-30  hours;  this  corre¬ 
sponds  to  the  period  of  transient  heaving. 

A  discussion  about  pore  water  migration  during  transient  heaving,  the  deformation 
activity  shown  in  Figure  17,  and  the  temperature-dependent  unfrozen  water  content 
listed  in  Table  2  has  been  published  by  Akagawa  (1988a)  regarding  the  location  and 
dominance  of  in-situ  and  segregation  freezing. 

Apparent  temperature  range  of  the  frozen  fringe 

From  the  temperature  profiles  and  the  location  of  the  warm  side  of  the  warmest  ice  lens 
shown  in  Appendices  C  and  E,  the  apparen  t  tempera  ture  of  the  segregating  ice  lens  Tg  was 
determined  (Fig.  18).  During  Test  A,  tended  to  rise  with  time,  from  -1°  to  -C.4°C,  during 
constant  overburden  pressure  and  boundary  temperatures;  during  Test  B,  Tg  tended  to 
remain  constant  ai  slightly  below  -0.8°C  during  constant  test  conditions.  Figure  19  shows 
the  difference  in  T$  between  Test  A  and  Test  B  when  the  test  conditions  were  kept  constant. 
It  is  uncertain  whether  there  is  a  difference  in  T$  in  the  first  30  hours  because  the  accuracy 
of  the  T$  values  for  the  the  warm  side  of  the  warmest  ice  lens  is  questionable. 

After  379  hours  in~est  A  and  383  hours  in  Test  B  the  overburden  pressure  was  changed 
several  times,  and  after  360, 377, 930.5  and  1006  hours  in  Test  A  and  1 1 1  and  1 67  hours  in 
Test  B  the  boundary  temperatures  were  changed  to  see  the  effect  on  T .  However,  the  data 
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Figure  17.  Strain  rate  development  in  the  freezing  soil.  Heave  is  positive. 


15  (bottom) 


(-•M/%) 
ajoy  uidjjs 


26 


Figure  1 7  (cont'd). 


a.  Test  A. 


b.  Test  B. 
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Figure  19.  Comparison  of  T.  in  Test  A  and  Test  B. 


points  in  this  report  might  not  be  adequate  to  describe  a  general  trend  for  the  effects  of 
overburden  pressure  and  boundary  temperature  on  T&.  Therefore,  the  data  are  shown 
only  for  reference. 

Apparent  thickness  of  the  frozen  fringe 

From  the  0°C  isotherm  depth  and  the  location  of  the  warmest  ice  lens,  shown  in 
Appendices  C  and  E,  the  apparent  thickness  of  the  frozen  fringe  L((  was  observed  (Fig.  20). 
During  Test  A,  Lff  tended  to  become  thinner  with  time,  from  1 .1  to  0.5  cm,  during  constant 
overburden  pressure  and  boundary  temperatures;  during  Test  B,  Lff  showed  the  same 
tendency  as  in  Test  A  for  the  first  30  hours  but  thenseemed  to  stay  at  about  1  cm  during 
constant  test  conditions  between  30  and  100  hours.  Figure  21  shows  the  difference  in  L(( 
between  Test  A  and  Test  B.  Again,  the  data  may  not  be  adequate  to  describe  the  effects  of 
overburden  pressure  and  boundary  temperature  on  L({. 

Thermal  conductivity  of  the  frozen  fringe 
during  stationary  heaving 

During  stationary  heaving  the  latent  and  sensible  heat  at  the  0°C  isotherm  are 
negligible.  The  thermal  conductivity  of  the  frozen  fringe  during  stationary  heaving  can  be 
demonstrated  as  follows.  The  heat  balance  at  the  0°C  isotherm  where  no  ice  growth  is 
expected  is 

Q«  =  Off  (3) 

where  Qu  and  Qff  are  heat  fluxes  in  unfrozen  soil  and  the  frozen  fringe,  respectively. 
Equation  3  can  be  rewritten  in  terms  of  thermal  conductivity  and  temperature  gradient  as 
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Figure  21.  Comparison  of  in  Test  A  and  Test  8. 
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Ku(d0/dz)u  =Kff(d0/dz)/f  (4) 

where  K  is  the  thermal  conductivity  of  thesoil  and  (d0/dz)  is  the  temperature  gradient  of 
thesoil.Thesubscripts  "u"  and  "ff"  stand  for  unfrozen  soil  and  frozen  fringe,  respectively. 
The  thermal  conductivity  of  the  unfrozen  soil  has  been  measured  by  the  thermal  probe 
method  (Penner  1975),  as  shown  in  Table  2.  The  temperature  gradients  within  the 
unfrozen  soil  (d0/dz)u  and  the  frozen  fringe  (d0/dz)ff  were  calculated  from  the  data 
listed  in  Appendix  A.  The  heat  flow  density  KJidQ/ dz)u  in  unfrozen  soil  was  determined 
by  using  the  measured  thermal  conductivity  from  Table  2  and  the  mean  temperature 
gradient  in  unfrozen  soil  calculated  from  Appendix  A. 

Most  of  the  temperature  gradients  in  the  frozen  fringe  were  calculated  using  the 
coordinates  and  temperatures  of  two  lead  spheres.  These  spheres  were  located  below  the 
0°C  isotherm  and  above  the  warmest  ice  lens  in  a  zone  where  no  expansive  strain  was 
observed.  However,  since  the  frozen  fringe  was  too  thin  to  contain  the  two  spheres 
between  44.5  and  832.5  hours  in  Test  A,  (d0/dz)ff  was  calculated  by  another  procedure 
for  this  layer.  First  the  0°C  isotherm  depth  was  extrapolated  by  assuming  a  linear  tempera¬ 
ture  profile  within  the  unfrozen  soil.  By  using  the  depth  of  the  0°C  isotherm  and  the 
temperature  and  coordinates  of  a  lead  sphere  located  in  the  shrinking  zone  in  the  frozen 
fringe,  (d0/dz)ff  was  calculated  for  this  period.  In  this  way,  the  thermal  conductivity  of 
the  frozen  fringe  K((  where  there  was  no  expansive  strain  was  calculated  by  substituting 
values  for  Ku,  (d0/dz)u,  and  (d0/dz)ff  in  eq  4: 


Kff  =Ku[(d0/rfz)u/(d0/dz)ff].  (5) 

Again,  since  the  0°C  isotherm  was  no  longer  penetrating  into  the  unfrozen  soil  (i.e.  there 
was  no  temperature  change  at  the  vicinity  of  the  0°C  isotherm)  in  this  stationary  heaving 
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Time  (hr) 


a.  Test  A. 


Time  (hr) 


b.  Test  B. 


Figure  22.  Heal  flow  density  in  the  unfrozen  soil  and  thermal  conductivity  of  the  frozen  fringe.  K  ( is  the  thermal 
conductivity  of  the  unfrozen  soil  observed  by  the  thermal  probe  method. 
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phase  and  no  ice  growth  in  the  shrinking  zone  was  expected,  the  latent  heat  and  the 
sensible  heat  at  the  0°C  isotherm  were  neglected  in  this  analysis. 

Figure  22  shows  the  calculated  heat  flow  density  and  thermal  conductivity  in  unfrozen 
soil.  The  heat  flow  density  in  unfrozen  soil  increases  with  time  up  to  about  30  hours  and 
then  stabilizes.  After  the  first  30  hours,  during  a  period  when  the  overburden  pressures 
were  maintained  at  0.06  MPa  for  Test  A  and  0.1 1  MPa  for  Test  B,  the  thermal  conductivity 
of  the  frozen  fringe  in  Test  A  decreased  greatly  from  the  value  for  unfrozen  soil  (KJ, 
whereas  there  was  little  decrease  from  Ku  in  Test  B.  However,  the  expected  error  in 
temperature  measurement,  shown  as  the  length  of  the  bar  at  each  data  point,  suggests  that 
more  data  are  needed  for  a  general  understanding  of  the  thermal  properties  during 
freezing. 

After 329  hours  in  Test  Aand  111  hours  inTest  B,  the  overburden  pressure  and  thecold- 
side  temperature  were  changed  several  times.  However,  there  was  no  clear  correlation 
with  the  thermal  conductivity. 

In  conclusion,  the  thermal  conductivity  in  the  shrinking  zone  of  the  frozen  fringe 
during  the  growth  of  the  final  ice  lens  is  similar  to  or  lower  than  that  of  the  unfrozen  soil 
(Fig.  22). 

Problems  in  determining  Ts  and  L({ 

So  far,  T  has  been  assumed  to  be  the  temperature  at  the  warm  side  of  the  visible 
warmest  ice  lens.  However,  much  of  the  data  in  Appendix  E  may  not  support  this 
assumption,  especially  during  transient  heaving.  Many  figures  in  Appendix  E,  which 
demonstrate  strain  rate  distributions,  show  the  expanding  layer(s)  between  the  locations 
of  the  warmest  ice  lens  and  the  0°C  isotherm,  especially  during  transient  heaving.  The 
existence  of  this  expanding  layer  may  mean  that  the  visible  ice  lens  does  not  indicate  the 
real  ice  segregation  location  but  only  the  result  of  the  ice  segregation. 

Another  empirical  result  that  supports  this  understanding  is  shown  in  Figure  23 
(Akagawa  1988b),  obtained  from  the  intensity  profile  of  the  X-ray  photos.  This  figure 


50  40 

Depth  from  the  Warm  Plate  (cm) 

a.  Test  A. 


Figure  23.  Difference  in  frozen  fringe  location  determined  from  the  location  of  the  warm 
side  of  the  warmest  icelensand  from  the  expansion  rate  profile  converted  from  the  intensity 
profile  of  the  X-ray  photo. 
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Consolidation  Heave 


60  50  40  30 

Depth  from  the  Warm  Plate  (cm) 

c.  Test  C. 

Figure  23  (coni' d). 

shows  the  continuous  strain  field  over  the  specimen,  along  with  the  locations  of  the  0°C 
isotherm  and  the  warm  side  of  the  visible  ice  lens.  The  location  of  the  warm  side  of  the 
visible  ice  lens  is  somewhat  behind  the  location  where  expansion  is  taking  place. 
Unfortunately  it  is  not  certain  what  is  causing  the  expansion  in  this  area.  The  expansion 
may  be  caused  by  the  early  stages  of  the  segregation  freezing,  but  there  is  no  proof  of  this. 

By  analyzing  the  shape  of  the  intensity  profile  between  the  0°C'  isotherm  and  the  location 
of  the  warmest  ice  lens,  it  might  be  possible  to  define  what  is  causing  the  expansion. 


33 


However,  until  it  becomes  clear,  the  determination  of  the  T$  and  Lffhas  to  rely  on  the 
location  of  the  warm  side  of  the  visible  warmest  ice  lens,  and  Tg and  L({ have  to  be  called 
"apparent  Tgand  L(f."  This  is  the  reason  that  qualitative  or  theoretical  analyses  have  not 
been  presented  in  this  report. 


SUMMARY  AND  CONCLUSIONS 

The  X-ray  radiograph  technique  is  effective  for  studying  freezing  soil  behavior, 
especially  when  the  lead  spheres  contain  thermocouples.  The  lead  spheres  enable  us  to 
determine  the  temperature  profiles,  and  the  X-ray  photographs  clearly  show  the  ice  lens 
distribution.  This  information  allows  us  to  determine  the  apparent  segregation  tempera¬ 
ture  and  the  thickness  of  the  frozen  fringe.  From  the  sphere  coordinates,  the  strain  and 
strain  rate  for  each  soil  layer  be  tween  ad  jacen  t  lead  spheres  can  be  calculated,  and  then  the 
distributions  of  strain  and  strain  rate  over  the  specimen  can  be  determined. 

The  tests  showed  that  the  frost  heave  activity  in  the  soil  while  freezing  is  as  follows. 
Before  the  0°C  isotherm  penetrates  through  a  soil  layer,  the  soil  layer  shrinks  and  then 
expands  during  transient  heaving.  The  amount  and  duration  of  this  shrinkage  and 
expansion  seem  to  be  closely  related  to  the  rate  of  the  0°C  isotherm  penetration,  which  in 
turn  has  a  close  relation  with  the  net  heat  flow  at  the  freezing  front.  However,  simple 
comparisons  between  the  0°C  isotherm  penetration  and  the  total  heave  or  total  heave  rate 
may  cause  confusing  results  because  of  the  shrinkage.  During  this  period  when  the  0°C 
isotherm  is  penetrating  into  unfrozen  soil,  not  all  of  the  net  heat  flow  is  expended  for  ice 
lens  growth  (segregation  freezing);  some  is  expended  for  in-situ  freezing,  which  will  not 
result  in  ice  lens  formation.  (Evaluations  of  expansion  amounts  due  to  in-situ  and 
segregation  freezing  are  important  for  determining  the  correlation  between  the  expansion 
rate  and  the  0°C  isotherm  penetration  rate). 

The  maximum  amount  of  the  shrinkage  is  as  high  as  10%  in  strain,  well  beyond  the 
preconsolidation  strain  of  this  soil.  From  an  analysis  of  the  magnitude  and  distribution  of 
the  shrinkage  with  load  and  porosity  data  obtained  from  a  consolidation  test,  the 
magnitude  of  the  local  pore  water  pressure  depression  was  found  to  be  greater  than  0.4 
MPa.  This  magnitude  should  have  a  close  relation  with  the  frost  heaving  force.  Unfrozen 
soil  also  shrinks  during  transient  freezing;  this  might  be  caused  by  a  pore  water  pressure 
gradient  created  by  the  steep  pore  water  drop  at  the  segregating  ice. 

This  study  shows  thefollowing  limitations  for  studying  frost  heave  using  conventional 
frost  heave  tests.  Since  during  a  period  when  a  CPC  isotherm  is  penetrating  (transient 
heaving)  the  frozen  fringe  and  the  unfrozen  soil  may  shrink,  the  amount  of  shrinkage  has 
to  be  determined  in  order  to  describe  the  frost  heave  properties.  Some  lateral  heat  flow 
may  not  be  avoidable  during  the  test,  so  the  tempeiature  measured  at  the  side  of  the 
specimen  may  contain  errors  caused  by  radial  heat  flow.  Since  the  temperature  profile  in 
frozen  soil  is  not  linear,  the  assumption  of  a  linear  temperature  profile  may  result  in  an 
error  of  some  tenths  of  a  degree  in  the  temperature  measurements. 

The  test  has  the  following  limitations  for  the  experimental  study  of  the  frozen  fringe. 
The  direct  strain  measurement  using  position  markers  like  lead  spheres,  as  in  this  study, 
might  not  be  precise  enough  for  studying  the  frozen  fringe.  However,  if  the  conversion 
from  intensity  of  the  X-ray  photos  to  strain  proposed  by  Akagawa  (1 988b)  is  introduced, 
strain  observations  in  a  small  area,  such  as  the  frozen  fringe,  may  become  accurate 
enough  Yfter  the  structure  of  the  frozen  fringe  is  clearly  observed,  the  accuracy  of  the 
temperature  measurements  demonstrated  in  this  study  may  be  sufficient.  However, 
without  the  experimentally  proven  structure  of  the  frozen  fringe,  the  absolute  accuracy 
of  0.1  °C  seems  to  be  inadequate.  Therefore,  the  most  crucial  breakthrough  in  the  study  of 
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Position  1 


APPENDIX  C.  X-RAY  PHOTOGRAPHS.  This  appendix  includes  all  the  radiographs  taken  during  Test  A  and  Test  B.  The  radiographs  wer 


n  in  pairs,  one  each  from  Position  1  and  Position  2. 
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frost  heaving  will  be  the  study  of  the  frozen  fringe  structure:  the  distribution  of  water  in 
this  zone  and  its  state  relative  to  the  thermal,  mechanical  and  chemical  conditions. 
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APPENDIX  A:  Z  COORDINATES  AND  TEMPERATURES, 
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(mm) 

(°C> 

(mm) 

(°C) 

(mm) 

(°C) 

Tw 

84.94 

2.3 

85.25 

2.2 

85.47 

2.3 

85.71 

2.2 

86.97 

2.2 

1 

79.39 

1.8 

79.74 

1.7 

80.02 

1.8 

80.31 

1.7 

80.55 

1.7 

2 

76.45 

1.6 

76.81 

1.5 

77.08 

1.5 

77.38 

1.5 

77.66 

1.4 

3 

72.3 

1.2 

72.69 

1.1 

72.99 

1.1 

73.25 

1 . 1 

73.41 

1 

4 

66.09 

0.7 

66.45 

0.6 

66.72 

0.6 

67.01 

0.6 

67.26 

0.5 

5 

63.3 

0.4 

63.64 

0.2 

63.91 

0.3 

64.17 

0.2 

64.42 

0.1 

6 

58.25 

0 

58.56 

-0.1 

58.85 

-0.1 

59.13 

-0.1 

59.4 

-0.3 

7 

52.68 

-0.5 

53.02 

-0.6 

53.33 

-0.6 

53.61 

-0.7 

53.82 

-0.8 

8 

49.01 

-0.8 

49.07 

-0.9 

49.13 

-0.9 

49.23 

-1 

49.31 

-1.2 

9 

43.7 

-1.3 

43.73 

-1.4 

43.74 

-1.4 

43.78 

-1.5 

43.82 

-1.6 

10 

37.74 

-1.8 

37.78 

-2 

37.85 

-i .  9 

37.86 

-2 

37.91 

-2.2 

11 

32.8 

-2.3 

32.79 

-2.4 

32.84 

-2.4 

32.89 

-2.5 

32.91 

-2.6 

12 

28.42 

-2.9 

28.41 

-3 

28.4 

-2.9 

28.46 

-3 

28.6 

-3.1 

13 

22.25 

-3.4 

22.41 

-3.6 

22.5 

-3.5 

22.53 

-3.6 

22.63 

-3.7 

14 

11.9 

-4.6 

11.9 

-4.6 

11.9 

-4.6 

11.9 

-4.6 

11.9 

-4.7 

Tc 

0 

-5.8 

0 

-5.9 

0 

-5.8 

0 

-5.9 

0 

-5.9 

37 


TEST 


A  (CONT'D) 


Time 

(hr)  :  15. 

5 

16. 

5 

17. 

5 

18. 

5 

19. 

5 

Z 

Tmp. 

Z 

Tmp. 

Z 

Tmp. 

Z 

Tmp. 

Z 

Tmp. 

(mm) 

<°C) 

(mm) 

<°C) 

(mm) 

(°C) 

(mm) 

(°C) 

(mm) 

<°C) 

No. 

Tw 

86.22 

2.2 

86.34 

2.2 

86.49 

2.2 

86.72 

2.2 

86.98 

2.2 

1 

80.84 

1.7 

81.05 

1.7 

81.26 

1.7 

81.48 

1.7 

81.69 

1.6 

2 

77.91 

1.5 

78.03 

1.4 

78.18 

1.4 

78.42 

1.4 

78.69 

1.4 

3 

73.7 

1 

73.96 

1 

74.15 

0.9 

74.36 

0.9 

74.55 

0.9 

4 

67.55 

0.5 

67.79 

0.5 

68.05 

0.4 

68.31 

0.4 

68.51 

0.3 

5 

64.68 

0.1 

64.89 

0.1 

65.11 

0 

65.36 

0 

65.6 

-0.1 

6 

59.69 

-0.2 

59.93 

-0.3 

60.13 

-0.4 

60.36 

-0.4 

60.57 

-0.4 

7 

53.94 

-0.7 

54 

-0.8 

54.06 

-0.9 

54.11 

-0.9 

54.12 

-1 

8 

49.36 

-1 . 1 

49.36 

-1.2 

49.4 

-1.3 

49.44 

-1.3 

49.45 

-1.4 

9 

43.87 

-1.6 

43.87 

-1.7 

43.85 

-1.8 

43.88 

-1.8 

43.9 

-1.8 

10 

37.9 

-2.1 

37.92 

-2.2 

37.93 

-2.3 

37.99 

-2.3 

37.99 

-2.4 

11 

32.9 

-2.6 

32.79 

-2.6 

32.8 

-2.7 

32.87 

-2.7 

32.93 

-2.8 

12 

28.52 

-3.1 

28.51 

-3.1 

28.47 

-3.2 

28.48 

-3.2 

28.46 

-3.3 

13 

22.59 

-3.6 

22.59 

-3.7 

22.6 

-3.7 

22.63 

-3.8 

22.66 

-3.8 

14 

11.9 

-4.6 

11.9 

-4.7 

11.9 

-4.7 

11.9 

-4.8 

11.9 

-4.8 

Tc 

0 

-5.8 

0 

-5.9 

0 

-5.9 

0 

-5.9 

0 

-5.9 

Time 

(hr):  20. 

5 

21. 

5 

22. 

5 

23. 

5 

44. 

5 

Z 

Tmp. 

Z 

Tmp. 

Z 

Tmp. 

Z 

Tmp. 

Z 

Tmp. 

(mm) 

(°C) 

(mm) 

<°C) 

(mm) 

<°C) 

(mm) 

(°C) 

(mm) 

(°C) 

No. 

Tw 

87.21 

2.2 

87.38 

2.2 

87.56 

2.2 

87.65 

2.2 

90.48 

2.2 

1 

81.9 

1.7 

82.07 

1.7 

82.26 

1.7 

82.37 

1.7 

85.15 

1.6 

2 

78.93 

1.4 

79.1 

1.4 

79.29 

1.4 

79.39 

1.4 

82.32 

1.3 

3 

74.74 

0.9 

74.93 

0.9 

75.13 

0.9 

75.24 

0.9 

78.2 

0.8 

4 

68.68 

0.4 

68.85 

0.4 

69.04 

0.4 

69.15 

0.4 

72.06 

0.2 

5 

65.79 

0 

65.99 

0 

66.18 

0 

66.31 

0 

69.18 

-0.2 

6 

60.8 

-0.4 

60.98 

-0.5 

61.19 

-0.4 

61.3 

-0.4 

62.39 

-0.8 

7 

54.09 

-1 

54.12 

-1.1 

54.16 

-1.1 

54.21 

-1 

54.41 

-1.3 

8 

49.42 

-1.4 

49.42 

-1.4 

49.44 

-1.4 

49.49 

-1.4 

49.64 

-1.6 

9 

43.91 

-1.8 

43.89 

-1.9 

43.91 

-1.8 

43.94 

-1.8 

44.09 

-2.1 

10 

37.9 

-2.3 

37.82 

-2.3 

37.87 

-2.3 

37.97 

-2.3 

38.09 

-2.4 

11 

32.89 

-2.7 

32.81 

-2.8 

32.81 

-2.7 

32.83 

-2.7 

32.97 

-2.8 

12 

28.48 

-3.2 

28.44 

-3.2 

28.62 

-3.2 

28.73 

-3.2 

28.48 

-3.3 

13 

22.52 

-3.8 

22.47 

-3.7 

22.5 

-3.7 

22.62 

-3.7 

22.64 

-3.7 

14 

11.9 

-4.7 

11.9 

-4.7 

11.9 

-4.7 

11.9 

-4.7 

11.9 

-4.6 

Tc 

0 

-5.9 

0 

-5.9 

0 

-5.8 

0 

-5.9 

0 

-6 

Time 

(hr) :  141 

.5 

304 

.5 

524 

.5 

664 

.5 

832 

.5 

Z 

Tmp. 

Z 

Tmp. 

Z 

Tmp. 

Z 

Trap 

Z 

Tmp. 

(mm) 

<°C) 

(mm) 

<°C) 

(mm) 

(°C) 

(mm) 

<°C) 

(mm) 

(°C) 

No. 

Tw 

95.89 

2.3 

101.03 

2.4 

103.39 

2.2 

103.96 

2.2 

106.52 

1.9 

1 

90.62 

1.8 

95.81 

1.8 

97.67 

1.7 

98.25 

1.7 

100.87 

1.3 

2 

87.88 

1 .  4 

93.17 

1.5 

95.06 

1.4 

95.66 

1.3 

98.35 

1 

3 

83.55 

0.9 

88.32 

1 

91.3 

0.9 

91.11 

0.9 

94.34 

0.7 

4 

77.36 

0.3 

82.4 

0.3 

85.28 

-0.3 

85.31 

0.3 

88.31 

0.1 

5 

74.34 

-0.2 

79.43 

-0.1 

82.43 

-0.1 

82.33 

-0.1 

85.34 

-0.3 

6 

62.59 

-1 

62.67 

-1  .  1 

62.94 

-1.2 

62.89 

-1.2 

62.99 

-1.3 

7 

54.12 

-1.5 

54.19 

-1.6 

— 

— 

— 

— 

— 

— 

8 

49.85 

-1.8 

— 

— 

— 

— 

— 

— 

— 

— 

9 

44.23 

-2.2 

44.31 

-2.3 

45.05 

-2.3 

45.06 

-2.3 

45.25 

-2.3 

10 

38.01 

-2.7 

38.1 

-2.8 

38.89 

-2.8 

38.73 

-2.8 

38.72 

-2.8 

11 

32.98 

-3.1 

33.34 

-3.1 

33.77 

-2.9 

33.58 

-2.9 

33.56 

-2.8 

12 

28.61 

-3.5 

28.44 

-3.5 

29.26 

-3.5 

29.11 

-3.4 

29.11 

-3.4 

13 

22.29 

-3.9 

22.47 

-3.9 

23.54 

-4.1 

22.93 

-4 

22.66 

-4.1 

14 

11.9 

-4.8 

11.9 

-4.7 

11.9 

-4.7 

11.9 

-4.7 

11.9 

-4.7 

Tc 

0 

-5.8 

0 

-5.8 

0 

-5.9 

0 

-5.9 

0 

-5.9 

38 


TEST  B 


Time  (hr) :  Ref  0 

2  Tmp.  Z  Tmp. 

(mm)  <°C)  (mm)  (°C) 

No. 

Tw  98.27  3.2  97.83  3.1 

0  92.34  3.1  91.9  3 

1  87.18  3.1  86.76  2.9 

2  82.72  3.1  82.31  2.8 

3  7/. 91  3  77.66  2.7 

4  72.89  3  72.69  2.7 

5  69.6  3  69.35  2.6 

6  63.3  2.9  63.21  2.5 

7  57.65  2.9  57.5  2.4 

8  52.46  2.9  52.24  2  5 

9  49.29  2.8  49.19  2.1 

10  42.62  2.7  42.42  2 

11  38.23  2.7  38.13  1.9 

12  33.09  2.6  33.05  1.6 

13  27.76  2.5  27.71  1.2 

14  17-51  2.2  17.59  0.3 

Tc  0  2.2  0  0 


1  2 
Z  Tmp .  Z 

(mm)  (°C)  (mm) 

98.35  3  98.32 

92.43  3.1  92.4 

87.47  2.9  87.34 

83  2.8  82.95 

78.38  2.7  78.34 

73.25  2.5  73.36 

70  2.4  70.07 

63.93  2.2  64.1 

58.1  1.9  58.35 

52.82  1.9  53.23 

49.82  1.5  50.22 

43.04  1.3  43.85 

38.57  1.1  39 . 5 1 

33.51  0.8  34.62 

28.1  0.4  29.4 

18.02  -0.3  19.31 

0  -5.4  0 


3 

Tmp.  Z  Tmp. 

(°C)  (mm)  (°C) 

3.1  98.5  3 

3  92.57  2.9 

2.8  87.52  2.6 

2.6  83.15  2.4 

2.4  78.48  2.1 

2.2  73.56  1.9 

2.1  70.25  1.7 

1.8  64.25  1.4 

1.5  58.67  1.1 

1.2  53.51  0.8 

1.1  50.66  0.6 

0.9  44.31  0.4 

0.6  40.02  0.1 

0.3  35.12  -0.2 

-0.1  30.05  -0. 

-1  19.28  -1. 

-5.4  0  -5.4 


Time  (hr) :  4  5 


l 

7 

Tmp. 

Z 

Tmp. 

Z 

Tmp. 

Z 

Tmp. 

2 

Tmp. 

No. 

( mm » 

(3C) 

(mm) 

(°C) 

(mm) 

<°C) 

(mm) 

(°C) 

(mm) 

(°C) 

Tw 

98 

.58 

3 

99 

.04 

3 

99.25 

3.1 

99 

.53 

3.1 

99 

.85 

3 

0 

92 

.  65 

2.9 

93 

.12 

2.9 

93.33 

2.9 

93 

.6 

2.8 

93 

.92 

2.7 

1 

87 

.68 

2.6 

87 

.95 

2.5 

88.16 

2.6 

88 

.5 

2.4 

88 

.  78 

2 . 3 

2 

83 

.34 

2.3 

83 

.61 

2.2 

83.9 

2.3 

84 

.18 

2.1 

84 

.39 

2 

3 

78 

.  7 

2 

78 

.92 

1.9 

79.34 

2 

79 

.54 

1.8 

79 

.8 

1 . 7 

4 

73 

.86 

1  .8 

74 

.04 

1.7 

74.24 

1.7 

74 

.64 

1.5 

74 

.94 

1  .  4 

5 

70 

.  56 

1.5 

70 

.88 

1.5 

71.11 

1.5 

71 

.39 

1.3 

71 

.72 

1 . 2 

6 

64 

.  58 

1.2 

64 

.91 

1 

65.03 

1 

65 

.38 

0.9 

65 

.75 

0.8 

7 

58 

.99 

0.8 

59 

.34 

0.7 

59.67 

0.7 

59 

.96 

0.5 

60 

.29 

0.4 

8 

53 

.92 

0.6 

54 

.31 

0.4 

54.69 

0.4 

54 

.95 

0.2 

55 

.31 

0.2 

9 

51 

.06 

0.3 

51 

.47 

0.2 

51.82 

0.2 

52 

.05 

0 

52 

.47 

-0.2 

10 

44 

.74 

0.  1 

45 

.18 

-0.1 

45.6 

-0.2 

45 

.98 

-0.4 

46 

.45 

-0.6 

11 

40 

.41 

-0.1 

41 

.03 

-0.3 

41.35 

-0.4 

41 

.59 

-0.8 

41 

.87 

-0.9 

12 

35 

.67 

-0.4 

36 

.06 

-0.7 

36.22 

-0.9 

36 

.24 

-1.3 

36 

.  44 

-1.5 

13 

30 

.36 

-0.9 

30 

.43 

-1.3 

30.4 

-1.6 

30 

.41 

-1.9 

30 

.52 

-2  .  1 

1 4 

19 

.49 

-2.4 

19 

.38 

-2.7 

19.5 

-2.9 

19 

.29 

-3.1 

19 

.3 

-3.3 

Tc 

0 

-5.4 

0 

-5.4 

0 

-5.4 

0 

-5.4 

0 

-5.4 

Time 

(hr)  : 

9 

10 

11 

12 

13 

7 

Ci 

Tmp. 

Z 

Tmp. 

Z 

Tmp. 

Z 

Tmp. 

Z 

Tmp. 

No. 

(mm) 

<°C) 

(mm) 

(°C) 

(ram) 

(°C) 

(mm) 

(°C) 

(mm) 

(°C) 

Tw 

99. 

95 

3 

100.28 

3 

100.56 

3 

100.79 

3 

101 

.03 

3 

0 

94 . 

02 

2.7 

94. 

36 

2.8 

94.64 

2.8 

94. 

87 

2.7 

95. 

1 

2.6 

1 

88. 

98 

2.3 

89. 

21 

2.4 

89.38 

2.4 

89. 

68 

2.3 

89. 

88 

2.3 

2 

84 . 

65 

2.1 

84. 

87 

2.1 

84.98 

2.1 

85. 

41 

2 

85. 

59 

1.9 

3 

80. 

07 

1.7 

80. 

26 

1.7 

80.48 

1.7 

80. 

79 

1.6 

81  . 

05 

1.6 

4 

75. 

23 

1.4 

75. 

47 

1.4 

75.62 

1.3 

75. 

96 

1.3 

76. 

17 

1.2 

5 

71  . 

99 

1.2 

72. 

18 

1 .  1 

72.3 

1.1 

72. 

7 

1 

72  . 

88 

1 

6 

66. 

05 

0.7 

66. 

33 

0.7 

66.62 

0.6 

66. 

77 

0.6 

67. 

19 

0.5 

7 

60. 

65 

0.4 

60. 

93 

0.3 

61.13 

0.2 

61  . 

46 

0.2 

61  . 

74 

0.  1 

8 

55 . 

64 

0.1 

55. 

97 

-0.1 

56.16 

-0.  1 

56. 

5 

-0.2 

56. 

85 

-0.3 

9 

52  . 

76 

-0.2 

53. 

07 

-0.4 

53.35 

-0.4 

53. 

59 

-0.5 

53. 

97 

-0.6 

10 

46  . 

89 

-0.7 

47. 

14 

-0.7 

47.22 

-0.8 

47. 

43 

-0.9 

47. 

53 

-1.1 

1 1 

42. 

06 

-1.1 

42. 

05 

-1.1 

42.14 

-1.2 

42. 

19 

-1.3 

42  . 

25 

-1.5 

1 2 

36. 

41 

-1.6 

36. 

44 

-1.7 

36.5 

-1.8 

36. 

52 

-1.9 

36. 

64 

-2 

13 

30. 

49 

-2.2 

30. 

41 

-2.2 

30.5 

-2.3 

30. 

52 

-2.4 

30. 

59 

-2 . 5 

1 4 

19. 

44 

-3.3 

19. 

34 

-3.4 

19.35 

-3.4 

19. 

39 

-3.5 

19. 

46 

-3.6 

Tc 

0 

-5 . 4 

0 

-5.4 

0 

-5.4 

0 

-5.4 

0 

-5.4 

39 


tO  Oi 


TES 


B  (CONT'D) . 


Time 

(hr)  : 

14 

15 

Z 

Tmp. 

Z 

Tmp. 

(mm) 

<°C) 

(mm) 

(°C) 

No. 

Tw 

101.29  3 

101.39  2.9 

0 

95.36 

2.6 

95.47 

2.7 

1 

90. 18 

2.3 

90.24 

2.3 

2 

85.84 

1.9 

85.92 

2 

3 

81.22 

1.5 

81.54 

1.6 

4 

76.38 

1.2 

76.61 

1.2 

5 

73.11 

1 

73.32 

0.9 

6 

67.35 

0.5 

67.52 

0.4 

7 

61.93 

</.  1 

62.15 

0 

8 

57.04 

-0.3 

57.25 

-0.4 

9 

54.13 

-0.6 

54.38 

-0.7 

10 

47.6 

-1.1 

17.58 

-1  . 1 

11 

42.2 

-1.5 

42.21 

-1.5 

12 

36.54 

-2.1 

36.56 

-2  .  1 

13 

30.56 

-2.6 

30.6 

-2.6 

14 

19.4 

-3.6 

19.32 

-3.6 

Tc 

0 

-5.4 

0 

-5.4 

Time 

(hr)  :  19 

20 

Z 

Tmp. 

Z 

Tmp. 

(mm) 

(°C) 

(mm) 

<°C) 

No. 

Tw 

n 

102.2 

2.9 

102.28 

2 . 9 

i 

o 

96.27 

91.12 

2.6 

2.2 

96.35 
91 .09 

2.6 

2 . 1 

c 

o 

86.88 

1.8 

86.92 

1 . 8 

o 

A 

82.27 

1.5 

82.27 

1  .  4 

4 

77.49 

1  .  1 

77.52 

1 

5 

6 

7 

74.19 

68.44 

0.8 

0.3 

74.38 

68.56 

0.8 

0.3 

l 

o 

63.05 

-0.1 

63.19 

-0.2 

o 

q 

58.21 

-0.5 

58.35 

-0.6 

1  0 
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APPENDIX  D:  STRAIN  RATE  OBSERVED  IN  THE  SOIL  LAYERS. 


Positive  strain  rate  represents  heave  and  negative  strain  rate  represents  consolidation. 


0  10  20  30  40  50 


Time  (hr) 


55 


S t ro i n  Rale 


10 


J _ 1—0  -J 

30  40  50 


40  50 


12 


40  50 


Time  (hr) 


Time  (hr) 


56 


Strain  Rate 


Strain  Rate 


Strom  Rate  (%/hr)  Strain  Rqte  (%/hr)  Strain  Rate  (%/hr) 


APPENDIX  E:  STRAIN  RATE  DISTRIBUTION  IN  THE  SOIL  LAYERS. 
Positive  strain  rate  represents  heave  and  negative  strain  rate  represents  consolidation. 
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APPENDIX  F:  LISTINGS  OF  T  .  Lff  AND  LOCATION  OF  THE  0°C  ISOTHERM. 


CO 

H 


0CO-'-HCONOtOOCOCOC,5COC)O«OC,)f3OOCO(r5nOCCCOCO(DNOO«) 

£OCOCOCOOCOCOCOOCOCOOlCOCOCOOOCOCOOC>COcOO>P*OCOOCOCOt,'COP' 


00000000000000000000000000000000 

i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  !  i  i  i  i  i  i  i  i 


£ 

w*  o 
-J  ~ 


t^»-<iOOU)CMlOCOrcMCO- iCOC'JIOCVJID— <Cv)*r*<r-D4CDOO*-<CO«-MtO*-HCSlCNJ 

I  I  I  MNOOONOOOOO^OOClOOOCOOOlOMXlOOOlOiOOCn 

II  I  *  •  . . 

<-r*-«~H^O-«Oi-<  —  o-*o»-<oo»-<^ooooo 


si 


OlOU)COwOWC'lHWTTnSSO'*<0'^tMU)'-<ClCOC5CMCOfMHHCOCOCO 

III . . . 

ocov,,r^'^'-^',rifiinioiowifli/)ioioiDioii5ti5toooooo0(OOcoo 


o 

CO  "s 

-  E 

o 

O 


rh-^'OCONMONNCOOlOSN^'OOlOCOlONCOOlOCO'HCOlNCOONSOlOO 

CMOCONSCO(MfO^SNO>OONC,3COC,3'^lfllOOy5ir3NHCOHN'Vin,<f'^,«r(OS 


(MCOCO,?'lJ",7lOlOlOlOlOlO(000(DOOOOOtOOOOSNNNSSNSSSS 


o  ^ 
E  ~ 

r*  ^ 


o 

^CJCO-'ClDlDS-COO— < 


CO  O  CM  -V  0  N  CO  CM  0 
NCOrWONCOClO^NCOTfCOCDO-HOtsHCOCOCOO^ 
IHfHrH^,^^Hr-(c'JCM(MCMWCM,^NH^^HCMWC,3COU)«D 


w 

H 


0)-^«3«>^COS^O)^SCOCO^OinSCMU5CMTfTrcOCOO) 

HOiOlSOODCOOOcOCOOSCOCOCOCOCOCOCOSCOWTrcocO’C 


— «ooooooooo*-«ooooooooooooooo 

I  I  I  I  I  I  t  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 


^  E 

Wh  O 

-J  ^ 


HHrjr)lOCMNCO|sO)-C)HNlfiSOCO{M^U5COOOlONN 
I  I  |^H-OHOO)OOCiNOiOOOOOCiOiOiCOWiflifliOifl 

III . . . 

<0— <00*m0^^-s-h-^000000000 


.J 

S  E  i 

O  | 

NJ  W 


h-Nr-'OS'ioiovcirHHcaiaHflcoh-ocoococor’TO 
i  Troo-Hrwscoa>oo*HCOCO^^iflu>(OM''{MO>cosSH 

I . 

fonr'CTrTr'winwirtjoifliowifiwinioioo^ssNco 


o 

cn  ■'*' 

-  E 

o 

o  ~ 


aiHiflOcocoNrioo>nco-<coNONHTr(ocoocMp)iOHifianDN 

No-i(ncoo)MroNcoooococoroiniooiooir>ooinoiNNN 


nco,r'<r'<r'?ifliniflifljfl<otooo(O0(fl0(DOoiO(Oi'*MHCOcocD 


o  ^ 
£  *— 

H  ~ 


Loinminiotoioioioioioioioioin 

in  ir>  to  in  lo  to  u>  m  uo  w  •  . . 

. CHcvHvrinioheooiC-iNCO^ 


in  m  in  to  in 


r  r  •?  cm 

TT  O  CO  CD  CD 
—<  CD  ID  CD  CO 


CO 

ON  H 


69 


I  so . : 1  oca t i on  of  the  zero  isotherm  from  the  cold  plate. 

WIL  :location  of  warm  side  of  the  visible  warmest  ice  lens 
from  the  cold  plate. 

: temperature  at  the  Z  WIL. 
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